Anisotropic collective flow of a Lorentz gas 



N Borghini and C Gombeaud 

Fakultat fur Physik, Universitat Bielefeld, Postfach 100131, D-33501 Bielefeld, 
Germany 

Abstract. Analytical results for the anisotropic collective flow of a Lorentz gas of 
massless particles scattering on fixed centres are presented. 



A remarkable feature of nucleus-nucleus collisions is the anisotropy of the particle 
emission pattern in the plane transverse to the collision axis: the transverse momentum 
distribution of outgoing particles reads 

d^A^ 1 dN 
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with Vnipr) the "anisotropic flow" coefficients, (p the azimuth of transverse momentum p 
and the event-by-event varying reference angle for the nth flow harmonic. Hereafter 
we shall neglect fluctuations, and all $„ will coincide with the x-axis. 

The experiment- driven focus of theoretical studies in the recent years has been 
on anisotropic flow for matter close to equilibrium. Here, we want to investigate the 
opposite case when particles undergo very few rescatterings, so that their evolution can 
meaningfully be described by a kinetic equation of the Boltzmann type. We specifically 
aim at obtaining analytical results — similar to those derived in P — which allow us to 
clearly identify qualitative behaviours together with their possible origins. 

As a further simplification, we consider the anisotropic fiow of a "Lorentz gas" of 
massless particles diffusing on infinitely massive particles. This constitutes a regular yet 
much simpler limiting case for the scattering of light particles on massive ones [2]. 

We wish to stress that the qualitative features which we derive in the following 
Sections are to our eyes more robust and thereby more important than the quantitative 
results. The model of a Lorentz gas may have little relevance for the phenomenology 
of heavy-ion collisions, yet it allows us to exemplify how in a more realistic description 
one should naturally expect 

• the mixing of different flow harmonics; 

• the evolution of anisotropic flow in the absence of spatial asymmetry when some 
flow is already present; 

• the non-monotonic time evolution of anisotropic flow. 

Our simple model also shows that such complex qualitative behaviours are not the 
exclusive privilege of approaches assuming many rescatterings like (dissipative) fluid 
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dynamics, but can appear quite generally, and play a role either at very early times [2111] 
or around the kinetic freeze-out, as well as for the anisotropic flow of fragile states. 



1. Expansion of a Lorentz gas 

Consider a gas of massless particles, described by the distribution density /(t,r,p), 
that scatter elastically with the differential cross section ad on a distribution nc(r) of 
fixed scattering centres. We shall assume that the problem is two-dimensional, i.e. we 
focus on the transverse dynamics of the gas, so that ad has the dimension of a length. 
/ then obeys the Boltzmann-Lorentz kinetic equation 

dJif r, p) + V Vr/(t, r, p) = n,(r) c j dG ad(e) [/(t, r, p') - f{f r, p)], (2) 

with V the particle velocity and the scattering angle of the diffusing particle. 

Integrating Equation ([2]) over space, the gradient term disappears, and one finds 
the evolution equation for the particle momentum distribution d^A^/d^p. The latter can 
then be multiplied by cos rup, with ip the azimuth of p, and averaged over ip, yielding 
the evolution equation for the anisotropic flow harmonic f„. 

At vanishing cross section, the solutions to Equation ([2]) are the free-streaming 
solutions 

/(°)(t,r,p) = /W(0,r-vt,p), (3) 

which are entirely determined by the initial distribution at t = 0. 

In the following, we study small deviations / = f^'^^ + f^^^ with \f^^^\ ^ /'-^^ to 
these solutions — which corresponds to considering very few scatterings per particle — by 
injecting the free-streaming solution in the collision integral in Equation ([2]). Introducing 
the total elastic cross section 

a,i, = y dead (6) 
and the unintegrated kernel 

C(X,p)^ y"d2rn,(r)/W(0,r-X,p), (4) 



one finds 



dt 



ifv) 



y"dead(e)C(v't,p')-aei.C(vt, p) 



(5) 



d2p 

The scattering rate at time t is given by 

r(t) = y"d2pd2rn,(r)/(t,r,p)aei.c^aei.cy"d2pC(vt,p). (6) 

Integrated over time, this rate gives the total number of rescatterings ACcat., which for 
the consistency of our approach should be small. 

For the density of scattering centres and the density distribution of diffusing 
particles at the initial time t = 0, we assume Gaussian profiles in position space 

^-W = 7T-7r^exp , (7) 
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with Nc the total number of scattering centres, and 

/'°'(0...P)^^expf-^-j4). (8) 



Airm^Ry 2i?2 2RI 

where the initial momentum distribution /(p) is normalized to 27r, so that the integral 
of /(O, r, p) over space and momentum yields the total number of diffusing particles. 
For the sake of simplicity we consider identical radii R^, Ry for both distributions. Let 

+ ^^y-l-e' 



With the initial profiles ([7]) and (j8]), the unintegrated kernel (jlj) reads 

XHl + e) + Y^(l- e " 



(9) 



with X = 



2. Isotropic initial momentum distribution, isotropic cross section 

Let us first consider the simplest case of an isotropic initial momentum distribution 
/(p) = /o(pt) as well as an isotropic differential cross section a a- The latter can then 
be replaced by adj^n and taken out of the gain term in Equation ([5]). Note that our 
normalization choice for /(p) is equivalent to 

j'OO _ 

/ fo{pT)PT<^PT = 1- 

Jo 

Let ip (resp. ip') denote the azimuth of p (resp. p') with respect to the direction of 
the X-axis of the scattering centre distribution ric- Then 

C(v«.P) = ^^^0^VT^e-^"'--\.,(-ff,ecos2^) . (10) 

and an analogous equation for C(v't, p'). 

This expression is readily integrated over p, yielding the rate 

with Iq the modified Bessel function of the first kind. Integrating from t = to infinity 
gives the total number of rescatterings over the evolution: 



where K denotes the complete elliptic integral of the first kind. At given Nc-, N, R and 
(Tel., this number of rescatterings is maximal for e = 0: one can thus fix the average 
number of rescatterings per diffusing particle f/scat./N at some small value in central 
collisions — which amounts to fixing the ratio Nc(Je\./R — and ensure a small number of 
rescatterings over all centralities. 
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The time evolution of the nth anisotropic flow harmonic follows from 



dtVn{t,p) = dt 



Jo d</^d^(^'P 



27r 



dipdt 



(^,P) 



cos rup, (13) 
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where in writing the denominator we have used the fact that elastic collisions on flxed 
scattering centres leave the momentum modulus pt unchanged, so that the momentum 
spectrum is actually independent of time. The integrand in the rightmost expression 
can be rewritten using Equation ([5]) and corresponds to the collision integral. The 
unintegrated kernel involved is given by Equation f lTOj) . First, the gain term depends on 
ip', not on ip, and thus does not contribute to dtVn{t,pT)- Then, the loss term yields 

for odd n; 



dtVn(f,PT) 



loss 



Nc(Tel.C 
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for even n. 



(14) 



Note that a factor of 2 is missing in the denominator of the equation as written in [2]. 

All V4n+2 coefficients, in particular f2, are increasing with time, while the Fourier 
harmonics f4„ are decreasing. Since the coefficients vanish at t = — the momentum 
distribution / is isotropic — one deduces for instance V2 > 0, but t>4 < 0: this reflects 
the alternating signs of the corresponding moments (in position space) of the initial 
Gaussian proflles. 

Integrating Equation f|T4l) over time from to t yields the time dependence of the 
Fourier coefficients Vnit^pr)- At early times |f„(t,pr)| oc f^^^, while at late times one 
finds for even n [5l formula 2.15.3(2)] 

VniPr) = }imyn(t,pT) 

" 11 



t—>-oo 



liVeael.v/r^ 



[n 



2"(f) 



2F, 



n + 1 n + 3 n r, 

■ — We 

4 ' 4 '2^ ' 



(15) 



4yii? 

where (2A; — 1)!! = 1 ■ 3 ■ ■ ■ (2A; — 1) if /c > 1, 1 if A; = 0, while 2-F1 denotes the Gaussian 
hypergeometrical function. For n = 2 (resp. n = 4), this formula reduces to Equation 
(C4) (resp. (C5)) of [2J. Interestingly, Vnipr) scales as e"/^ for small eccentricities. 

In Figure [1], we show the dependence of V2 on impact parameter b in Pb-Pb 
collisions, where the eccentricity e is related to b through the Glauber optical model, 
assuming that the eccentricity dependence of V2 is given by Equation (fT5|) with n = 2. 



3. Anisotropic initial momentum distribution, isotropic cross section 

We now allow for the possibility that the expanding gas possess an initially anisotropic 
momentum distribution, which we describe by introducing its Fourier series 



/(p) = Mpt) 



1 + 2 {wk^c cos kip + Wk^s sin kip) 



k=l 



(16) 



The Fourier coefficients Wk^c, Wk,s could generally depend on p^, yet we shall hereafter 
leave this dependence aside. The coefficients Wn,c of the cosine harmonics correspond to 
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b (fm) 



Figure 1. Centrality dependence of V2 (here, for an evolution with on average 0.1 
coUision per particle). 



the "usual" anisotropic flow coefficients, taken at the initial time, Vn{t = 0,pt)- As in 
the previous Section, the differential cross section is taken to be isotropic. 

The initial momentum distribution ( |T6l) gives for the unintegrated kernel 

NcNfoipT) 



C(vt,p) 



■VT 



g2 Q-cH^/iR" 



X exp 



e cos 2(p 



2 ^ {wk,c cos kip + Wk,s sin k(p) 



k=l 



.(17) 



With this kernel, the scattering rate is given by 



m 



g2 ^-cH^/iR^ 



and the total number of rescatterings, which has to be kept small, by 



AC 



NrNa. 



cl. 



scat. 



(2g-l)!! 
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The unintegrated kernel (JTTll also allows one to compute the time derivative of the 
anisotropic flow coefficient f„. As in Section |2l the gain term of the collision integral 
does not contribute to dtVn(t,pT), whereas the contribution of the loss term follows from 
multiplying Equation (fT7|) with cos mp and then integrating over ip. The only terms from 
the sum over k that result in a non- vanishing integral are those in cos k(p with k of the 
same parity as n, so that n — k and n + k are even. The isotropic part of the momentum 
distribution only contributes when n is even, as in Equation ([1] 
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All in all, one finds 



dtV2mit,PT] 



loss 
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Let us shortly discuss these results, focusing first on the short-time behaviour. 
Taking m = 1, Equation f llSap gives for the evolution of elliptic flow 



dtV2(t,pT) 



for t < — . 

c 



That is, in the presence of a positive initial elliptic flow, f 2(t, pr) first decreases (linearly), 
before it starts increasing: since more particles are emitted in-plane than out-of-plane, 
there are more particles "lost" at = or 180° than at ±90°. On the contrary, a 
negative initial W2,c = 'f^2(0,Pr) accelerates the initial increase of elliptic flow, while the 
later behaviour depends on the sign of e(l + W4^c) + 2ti'2,c- In either case, V2 evolves even 
for vanishing eccentricity e = 0, which is quite a nontrivial finding. One also sees that 
V2 is influenced by the presence of any finite initial w^^^c = V/i{f),pT)- 

For odd harmonics ra. Equation (118^1) shows that some finite Vn{t,PT) can develop 
if and only if there exist at least one non- vanishing odd harmonic at t = 0. In the case 
of directed flow Vi for instance. Equation (\18b\i with m = gives 

Nc(Tcl.C 



dtVi{t,pT) 



4i?2' 



167ri?4 

whereas for triangular flow v^, one finds 



Vl^ wi,cet^ + O(t') for t < -, 



dtV3{t,PT) 



for t < — . 

c 



Thus, f3 evolves even in the absence of any "triangularity" in the collision geometry — 
in obvious similarity to the evolution of V2 for e = 0. Additionally, dtv^it.px) again 
illustrates the mixing of different harmonics present in Equations ( ll8al) -( fT86jl . 
The latter can be integrated from t = to oo. One in particular gets 



V2{Pt) = W2,c + 



+ 2. T 



^ 5 2 



W2. 
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The initial elliptic flow W2,c = V2{t = 0, pt) breaks the linear scaling of f 2 with eccentricity 
at small e both trivially as well as through its influence on the anisotropic flow developed 
in the rescatterings: 



Again, we find the mixing of different harmonics as well as an evolving f 2 at e = 0. 
Note that the ratio N(.(Te\./ ^\/t^R necessarily takes a small value when the mean number 
of rescatterings per particle is small. Accordingly, V2{pt) does not differ much from its 
initial value, which is normal within our few-rescatterings approach. 

4. Isotropic initial momentum distribution, anisotropic cross section 

We now come back to an isotropic initial momentum distribution, but consider the case 
of an anisotropic differential cross section (Td(0). The latter can generally be expanded 
as a Fourier series. The first harmonic in the expansion describes an asymmetry 
between forward and backward scattering, the former being favored if the corresponding 
coefficient is positive. Then, the second harmonic accounts for increased or suppressed 
scattering at ±90° with respect to or 180°. Higher harmonics describe less natural 
behaviours, which we shall not consider in the following. Additionally, we assume that 
the interaction preserves parity, so that sine harmonics vanish. We thus restrict ourselves 
to a differential cross section given by 



Note that the coefficients and <^2 are not totally arbitrary, since cxd must remain non- 
negative when spans the range [0, 2tt]: one for example easily checks that, irrespective 
of the value of ^1, one should have |<^2| < |. 

The anisotropy of the differential cross section does not affect the scattering rate 
nor the resulting total number of rescatterings, which are thus given by Equations (ITTll 
and f|T2|) . The loss term of the collision integral relies on the total elastic cross section 
and is thus the same as in Section [21 it still yields a contribution to dtVn{t,pT) given by 
the right-hand side of Equation ([14]) . 

On the other hand, the gain term of the collision integral now gives a non-vanishing 
contribution, since ip' is no longer arbitrary, but related to (p through = + 6, with 
a non-uniform distribution in 9. Inspecting Equations ([5]), f[TO|) and f [T3|) together with 
the differential cross section ( |T9l) . the contribution to dtVn(t,p) of the gain term reads 





(19) 




gain 
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Irrespective of the value of n, the term leads to a vanishing integral over (f', while 
the integrals of the constant and ^2 terms yield modified Bessel functions, so that the 
expression between curly brackets equals 



while it vanishes for n ^ 2, i.e. the gain term only contributes to the second harmonic 
of anisotropic flow, that is elliptic flow. Putting the gain and loss terms together, one 
eventually obtains after integrating over time 



while Vnipr) for even n 2 remains given by Equation f lTSj) . Thus, an increased 
(resp. decreased) scattering probabihty at ±90°, as found e.g. in collisions of identical 
bosons (resp. fermions) — which is obviously not the case of the colliding particles in our 
Lorentz-gas model — , leads to a larger (resp. smaller) V2- 

Eventually, one can mix the various ingredients together and consider an anisotropic 
differential cross section together with some initial anisotropic flow. In that case, the ^1 
coefficient starts playing a role when combined with non-vanishing initial Wn,c, while ^2 
will affect further flow harmonics besides V2- 
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In turn, the remaining integral over (p is trivial and yields for n = 2 






